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ABSTRACT 

Aims. We study dark matter content in early-type galaxies and investigate whether X-ray luminosities of early-type galaxies are 
determined by the surrounding gravitational potential. 

Methods. We derived gravitational mass profiles of 22 early-type galaxies observed with XMM-Newton and Chandra. 
Results. Sixteen galaxies show constant or decreasing radial temperature profiles, and their X-ray luminosities are consistent with 
kinematical energy input from stellar mass loss. The temperature profiles of the other 6 galaxies increase with radius, and their X-ray 
luminosities are significantly higher. The integrated mass-to-light ratio of each galaxy is constant at that of stars within 0.5 ~ \r e , and 
increases with radius, where r e is the effective radius of a galaxy. The scatter of the central mass-to-light ratio of galaxies was less in 
K-band light. At 3r e , the integrated mass-to-light ratios of galaxies with flat or decreasing temperature profiles are twice the value at 
0.5r e , where the stellar mass dominates, and at 6r e , these increase to three times the value at 0.5r e . 

Conclusions. This feature should reflect common dark and stellar mass distributions in early-type galaxies: Within 3r e , the mass of 
dark matter is similar to the stellar mass, while within 6r e , the former is larger than the latter by a factor of two. By contrast, X-ray 
luminous galaxies have higher gravitational mass in the outer regions than X-ray faint galaxies. We describe these X-ray luminous 
galaxies as the central objects of large potential structures; the presence or absence of this potential is the main source of the large 
scatter in the X-ray luminosity. 
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1. Introduction 

The bottom-up hierarchical theory of galaxy formation predicts 
that gala xies should be emb edded in massive dark matter ha- 
los (e.g. iNavarro et alj|1997l) . The presence of dark matter in 
spiral galaxies has been revealed through observations of stel - 
lar rotation curves dRubin et all Il978t Ivan Albada et all Il985l) . 
However, the study of halos in early-type galaxies is limited due 
to a lack of suitable and easy tracers such as rotation curves. 
Recently, observations of stellar velocity dispersion of early- 
type galaxies have reached to 1 ~ 2 r e , and a correlation 
betw een mass-to-light ratio and optical luminosit y was found 
(e.g., iKronawitter et al.ll2000t iGerhard et alj|200ll) . Here, r e is 
the effective radius of a galaxy. For a small number of galax- 
ies, mass profiles up to several r e have been obtained using test 
particles such as globular clusters or planetary nebulae (e.g . , 
iRomanowskv et alJl2003l;IChakrabartv & Ravchaudhury|2008l 

X-ray observations provide a powerful tool to study the 
shape of the gravitational potential, and hence dark matter dis- 
tributions, of early-type galaxies. These galaxies have a hot, X- 
ray-emitting interstellar mediu m (ISM), which is considered to 
be gravitationally confined ( e.g., lForman et al.lll985t iMatsushital 
l200UlFukazawa et al.l2006l) . The ISM luminosities of early-type 
galaxies vary by two orders of magnitude for the same optical B- 
band luminosity (Lb) (e.g., [C anizares et al. 198^ iBeuing et al.l 
119991: IMatsushita et al.l 120001: IMatsushital 1200 lb . whereas opti- 
cal observations indicate that these galaxies are dynamically 
uniform systems (iDiorgovski & D avis 1987; lBender et al.ll 1993t 



IKronawitter et al.l 12000). A key to solving this discrepancy is 
the extended X-ray emissions that hav e been detected around 
X-ray luminous early-type galaxies (IMatsushita et al.l 19981: 
IMatsushital 1200 lb . On the basis of ROSAT data, IMatsushital 
(1200 ll) discovered that, for most early-type galaxies, ISM lu- 
minosities within the optical radius agree with kinematical en- 
ergy input from the stellar mass loss L a . These galaxies have 
flat or decreasing temperature profiles against radius. By con- 
trast, galaxies with ISM luminosities much larger than L lT show 
largely extended emission with a radius of a few tens of r e 
with positive temperature gradients. XMM-Newton RGS obser- 
vations provided evidence of a weak positive temperature gra- 
dient in the inner region of the ISM i n NGC 4636, w hich has a 
much higher ISM luminosity than Lg- dXu et al. 1 12051) . The cor- 
relation between the temperature gradient and spatial distribu- 
tion w as confirmed with Chandra observations dFukazawa et al.l 
120061) . These features suggest that X-ray luminous early-type 
galaxies commonly sit in the center of a large-scale (a few hun- 
dred kpc) potential well, which leads to their high luminosities. 
Other galaxies may lack such a large-sca le potential. On the ba- 
sis of the extent of the ISM brightness, IMatsushital d2001l) de- 
noted galaxies as either X-ray extended galaxies or X-ray com- 
pact galaxies. The gravitational mass profile of cD galaxies also 
shows two distinct contributions that can be assigned to the 
gravitational potential of the cD galaxy and that of the cluster 
(IMatsushita et al.l2002l) . Thus, the only way to measure the grav- 
itational mass profile of pure early-type galaxies is to observe the 
X-ray compact galaxies. 



2 



R. Nagino and K. Matsushita: Gravitational potential and X-ray luminosities of early-type galaxies 



With ROSAT PSPC observations, |Q' Sullivan etaT] d2003l) 
found that the relation between the central stellar velocity disper- 
sion and the temperature obtained from X-ray emission is similar 
to that for clusters and the relation between the X-ray luminos- 
ity and the temperature has a steep slope comparable with that 
found for galaxy groups. 

Chandra and XMM-Newton have already observed several 
tens of early-type galaxies. Most of the analysis was done 
for X-ray luminous and extended objects, and the number of 
X-ray compact galaxies with accurately derived gravitational 
mass profiles is still limited. For X-ray l uminous galaxies, mass 
profiles are easily ob tained over I0r e dFukazawa et al.l 120061: 
iHumphrev et al.ll2006l) . Using XMM-Newton, even for several 
X-ray c ompact galaxies, mass profiles can be derived up to sev- 
eral r e dFukazawa et al.ll2006l) . and observed gravitational mass 
profiles of X-ray extended and X-ray compact galaxies are sim- 
ilar when plotted against radius in units of r2oo- T he dark matter 
profil e s are well described by the NFW model dNavarro et alj 
119961 119971) . which is based on numerical simulations as- 
suming cold dark matter (CDM) as well as galaxy clusters 
dFukazawa et alj|2006t IZappacosta et al.l 120061) . Chandra obser- 
vations suggested that the shape of the X-ray is ophotes is unre- 
lated to the shape of the gravitational potential dDiehl & Statlerl 
[20071 l2008h . 

In this study, we obtained gravitational mass profiles of 22 
early-type galaxies observed with XMM-Newton and Chandra 
to investigate whether X-ray luminosities of early-type galax- 
ies are determined by the surrounding gravitational potential and 
to study dark matter content in early-type g alaxies. Throughout 
this pa per, we adopt the solar abundances of lAnders & Grevessd 
d 19891) . Unless otherwise specified, errors are quoted at lcr con- 
fidence. 



2. Targets and Observations 

We analyzed archival data of 22 early-type galaxies with dis- 
tances less than 40 Mpc and B-band luminosities Lb > 10 101 L o 
observed with XMM-Newton. The v alues of Lg a nd the dis- 
tances to the galaxies are taken from iTullvl d!988h . The char- 
acteristics and observational log of the sample galaxies are sum- 
marized in Tables [TJ and [2] respectively. The sample includes 
15 elliptical and 7 SO galaxies. Eight are located in the Virgo 
Cluster, 2 are in the Fornax Cluster, and the others are either in 
the field or in small groups. All observations were carried out 
with MOS1, MOS2, and PN together. 

We also used Chandra data for 19 of the sample galaxies 
with good signal-to-noise ratios to derive the mass profile at their 
central regions. As summarized in Table|2] 15 galaxies were ob- 
served with ACIS-S, and 4 galaxies were observed with ACIS-I. 



3. Data Reduction 

3.1. XMM-Newton 

We analyzed MOS1, MOS2, and PN data of 21 galaxies. For 
NGC 4472, only MOS1 and MOS2 data were used, since PN 
data for this galaxy did not exist in the archive. We used 
XMMSAS version 7.0.0 for the data reduction. 

We selected events with FLAG = and pattern smaller than 
4 and 12 for the PN and MOS, respectively. A significant frac- 
tion of XMM-Newton observations is contaminated by soft pro- 
ton flares. To filter the flares, for each observation, we made a 
count rate histogram of each detector, fitted the histogram with a 
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Fig.l. Raw MOS (MOS1 + MOS2) spectrum at r=6-14' of 
NGC 4636 (black), and the adopted background spectrum (red). 



Gaussian, and selected times within 2.5cr of the mean of the his- 
togram. The total exposure times after screening the flare events 
are summarized in Table [2] 

The spectra were accumulated within rings centered on the 
center of each galaxy. Hereafter, we denote r as the projected 
radius from the galaxy center. We excluded point sources with 
MOS1 and MOS2 count rates larger than 0.01 count/sec. The 
edetect_chain command was used to detect point sources. The 
response matrix file and the auxiliary response file correspond- 
ing to each spectrum were calculated using SAS version 7.0.0. 

The background spectrum was calculated for each spec- 
trum by integrating blank sky data in the same detector regions. 
Among deep sky observations with the XMM, we selected data 
with the most similar background to that of each galaxy, after 
screening background flare events in the same way. Each back- 
ground agrees well with the data at higher energies, as shown in 
Figure[TJfor NGC 4636. 

Table|2]also summarizes total counts of MOS and PN within 
4r e centered on each galaxy. Here, an annular region, 10-14', 
from the center of each galaxy was used as a background after 
subtracting a blank sky data. The total counts of sum of those 
of MOS and PN have a wide range from 1000 to 600000. We 
analyzed projected annular spectra of all of the sample galaxies, 
while fittings of deprojected spectra were performed for thirteen 
galaxies with the total counts > 12000. 

3.2. Chandra 

Chandra data analysis was performed with the CIAO software 
package, version 3.3. We excluded time regions with a high 
background rate. We also eliminated point sources identified 
with the tool wavedetect. The outer region of each data set was 
subtracted as background. 



4. Spectral Analysis and Results 

4.1. Spectral fit 

4.1 .1 . Projected annular spectra 

To derive the gravitational mass profiles of individual galaxies, 
we need temperature and density profiles of the ISM. First, we 
fitted projected annular spectra centered on each galaxy from 
MOS (MOS1 + MOS2) and PN simultaneously, except NGC 
4472. To exclude possible emissions from our Galaxy and sur- 
rounding clusters, we also subtracted the spectrum in an annular 
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Galaxy 


Type" 


D" 


r c 

1 e 


log L B " 


a e 


N H 1 


Note 






(Mpc) 


(arcmin) 


(L B ) 


(km/s) 


(10 20 cnr 2 ) 




IC1459 


-5.0 


20.0 


0.58 


10.42 


311 


1.18 


X c 


NGC720 


-5.0 


20.3 


0.60 


10.34 


240 


1.54 


X c 


NGC1316 


-2.0 


16.9 


1.35 


10.78 


250 


1.89 


Xc, Fornax 


NGC1332 


-2.0 


17.7 


0.47 


10.22 


319 


2.23 


X c 


NGC1395 


-5.0 


20.0 


0.81 


10.33 


254 


1.99 


x E 


NGC1399 


-5.0 


16.9 


0.68 


10.31 


362 


1.34 


X E , Fornax 


NGC1549 


-5.0 


13.4 


0.78 


10.10 


213 


1.46 


X c 


NGC3585 


-5.0 


21.6 


0.60 


10.56 


227 


5.58 


Xc 


NGC3607 


-2.0 


19.9 


0.73 


10.41 


223 


1.48 


X c 


NGC3665 


-2.0 


32.4 


0.48 


10.54 


186 


2.06 


Xc 


NGC3923 


-5.0 


25.8 


0.83 


10.68 


241 


6.21 


Xc 


NGC4365 


-5.0 


16.8 


0.83 


10.40 


266 


1.62 


X c ,Virgo 


NGC4382 


-2.0 


16.8 


0.91 


10.64 


187 


2.52 


X c ,Virgo 


NGC4472 


-5.0 


16.8 


1.74 


10.92 


302 


1.66 


X E , Virgo 


NGC4477 


-2.0 


16.8 


0.63 


10.14 


175 


2.64 


X c ,Virgo 


NGC4526 


-2.0 


16.8 


0.74 


10.41 


260 


1.65 


X c ,Virgo 


NGC4552 


-5.0 


16.8 


0.49 


10.35 


262 


2.57 


X c ,Virgo 


NGC4636 


-5.0 


17.0 


1.48 


10.46 


208 


1.81 


X E , Virgo 


NGC4649 


-5.0 


16.8 


1.15 


10.74 


343 


2.20 


X c ,Virgo 


NGC5044 


-5.0 


38.9 


0.89 


10.60 


237 


4.93 


x E 


NGC5322 


-5.0 


31.6 


0.56 


10.86 


233 


1.81 


X c 


NGC5846 


-5.0 


28.5 


1.05 


10.66 


251 


4.26 


x E 



"Morphological type code from Tullv (1988). 
^Distance to the galaxy from TufJ^ (1988) . 

'Effective radius from RC3 Catalog (de Vaucouleurs et al.. 1991). 

"Total B-band luminosity from Tullv (1988) . 

'Central stellar velocity dispersion from Prugniel & Simien (1996). 

^Column density of the Galactic absorption from Dickey & Lockman (1990). 



region, 10'-14', from the annular spectra o f each galaxy. The 
fitting model is a sum of a vAPEC model (ISmith et all boOll) 
and a power-law model. The vAPEC model represents thin ther- 
mal emission from the ISM, and the power-law model repre- 
sents the contribution from unresolved low-mass X-ray binaries 
(LMXBs), where we fixed the power-law index at 1.6. Since 
Chandra observations found that total spectra of discrete sources 
in earl y-type galaxies are well described with this power-law 
model (iBlanton et all 1200 UlRandall et al.ll2004l) . The two com- 
ponents were subjected to a common absor ption with fixed col- 
umn d ensity, Nh, at the Galactic value from lDickev & Lockmanl 
(119901) . We organized heavy element abundances into three 
groups: the ff-element O group (O, Ne, and Mg), Si group (Si 
and S), and Fe group (Fe and Ni). The abundances of the three 
elemental groups were allowed to vary. For the innermost region 
of IC 1459, we add ed a power-law model from the central nuclei 
found by Chandra dFabbiano et alll2003l) . 

For brightest galaxies, NGC 4472, NGC 4636, NGC 4649, 
and NGC 5044, whose total counts within 4r e larger than 
200000, we performed spectral fitting on each annular region. 
In order to derive accurate temperature profiles of the ISM in 
the other galaxies with lower signal-to-noise ratio, the spectra 
of all annular regions were fitted simultaneously, where the ISM 
abundances were assumed to have common values. 

Table [3] and Figure A. 1 summarize the results of the spec- 
tral fitting, ISM temperature, abundance, and ISM luminosity. 
Figure [2] shows MOS and PN spectra of the innermost regions 
several representative galaxies. The spectra of X-ray faint galax- 
ies, whose 

total counts within Ar e are smaller than 12000 counts, pro- 
vide gradually smaller values of reduced-^ 2 with this single- 
temperature model (hereafter IT model) for the ISM. However, 



a representative spectrum of NGC 4382 whose ISM tempera- 
ture is ~ 0.4 keV, there are residual structures around 0.9 keV. 
The X-ray brighter galaxies show larger reduced-^ 2 . The spectra 
of galaxies with kT~ 0.6 keV, NGC 4636, NGC 720 and NGC 
3923, show common residuals at 0.7-0.9 keV (Figure|2]i. While, 
NGC 4649, which the galaxy with kT~ 0.8 keV, has different 
residual structures from these three galaxies. 

4.1.2. Deprojected Spectra 

To consider projection effects, we performed spectral fittings of 
deprojected spectra. For the data with high statistics, total counts 
of sum of those of MOS and PN within Ar e larger than 12000, 
deprojected spectra were calculated using "onion peeling" meth- 
ods by subtracting the contribution from the outer shell regions 
for all spectral componen ts, assuming the ISM is spherically 
symmetric as described bv lTakaha shi (2004). We then fitted the 
deprojected spectra with the IT model in the same way as in 
Section |4~TTT] 

Tableland Figure A.l summarize the results. 

The reduced x 2 reduced to ~ 1 except the innermost re- 
gions of three brightest galaxies, NGC 4472, NGC 4636, and 
NGC 4649. However, residual structures around 0.7-0.9 keV 
in the projected spectra still remain in the deprojected spectra 
fitted with the IT model (Figure [2]). We also fitted the spectra 
of innermost regions of the three brightest galaxies with a two- 
temperature vAPEC model for the ISM (hereafter 2T model). 
The results are summarized in Table[4] Then, reduced^ 2 reduced 
to 1.2-1.5, and even the 2T model gives similar residual struc- 
tures (Figure [2]). These discrepancies in the Fe-L energy range 
are also seen in the RGS spectrum of the X-ray luminous ellip- 
tical galaxy, NGC 4636 dXu et all |2002|) . Suzaku observations 
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Table 2. Observational log of the sample galaxies 



XMM-Newton Chandra 



Galaxy ObsID" exposure" total counts'" ObsID" ACIS exposure 

(ksec) (xlO 3 ) (ksec) 



IC1459 


0135980201 


25,26,21 


8.9,14 


2196 


ACIS-S 


54 


NGC720 


0112300101 


17,18,11 


5.6,7.1 


492 


ACIS-S 


38 


NGC1316 


0302780101 


51,62,29 


32,38 


2022 


ACIS-S 


26 


NGC1332 


0304190101 


54,54,41 


10,16 


4372 


ACIS-S 


54 


NGC1395 


0305930101 


43,46,31 


15,21 


799 


ACIS-I 


15 


NGC1399 


0012830101 


3,3,3 


7.1,10 


240 


ACIS-S 


43 


NGC1549 


0205090201 


8,8,6 


0.4,1.0 








NGC3585 


0071340201 


11,11,7 


0.6,0.9 


2078 


ACIS-S 


35 


NGC3607 


0099030101 


14,14,6 


2.5,1.9 


2073 


ACIS-I 


38 


NGC3665 


0052140201 


23,24,19 


2.3,4.2 


3222 


ACIS-I 


18 


NGC3923 


0027340101 


32,32,24 


11,18 


1563 


ACIS-S 


19 


NGC4365 


0205090101 


25,25,21 


3.8,7.1 


5923 


ACIS-S 


37 


NGC4382 


0201670101 


16,16,14 


3.0,5.6 


2016 


ACIS-S 


39 


NGC4472 


0200130101 


79,80 — 


300 — 


321 


ACIS-S 


37 


NGC4477 


0112552101 


13,13,7 


1.7,1.9 








NGC4526 


0205010201 


20,20,16 


1.7,3.0 


3925 


ACIS-S 


41 


NGC4552 


0141570101 


21,22,15 


10,15 


2072 


ACIS-S 


53 


NGC4636 


0111190701 


58,58,50 


210,320 


4415 


ACIS-I 


73 


NGC4649 


0021540201 


46,46,36 


84,120 


785 


ACIS-S 


31 


NGC5044 


0037950101 


17,17,8 


130,92 


3225 


ACIS-S 


82 


NGC5322 


0071340501 


16,15,12 


1.0,1.5 








NGC5846 


0021540501 


13,13,9 


25,28 


788 


ACIS-S 


24 



"Observation number of the XMM-Newton and Chandra data. 
'"Exposure time of the EPIC-MOS1, MOS2, and PN, respectively. 

c Total counts within 4r f at 0.3-2.0 keV of the MOS (MOS 1 + MOS2) and PM, respectively. 



of NGC 720 and NGC 1404 whose ISM temp eratures are also ~ 
0.6 k e V also give s i milar residual structures dMatsushita et al.L 
120071: iTawara et al.L 120081) . Therefore, these residual structures 
are likely to be related to poorly modeled Fe-L lines. 

The spectral fittings of the deprojected spectra give mostly 
same temperatures with those of the projected ones (Figure A. 1). 
Therefore, for the fainter galaxies which were not performed the 
deprojected analysis, we used the temperature profiles derived 
from the projected spectra to derive gravitational mass profiles. 

4.2. Results 

4.2.1. ISM luminosities within 4r c 

We derived the absorption-corrected ISM luminosities (Lism) 
and the luminosities of the power-law component (L/ M „/) in the 
energy band of 0.3-2.0 keV within Ar e (Table |5). The values 
of Lism derived from the deprojected spectra are close to those 
from projected annular spectra. Hereafter, we use Lism from the 
deprojected analysis. For the X-ray fainter galaxies without de- 
projected analysis, we use the values derived from the projected 
analysis. 

The relationship of Lism to Lb is shown in the left panel of 
Figure[3] The sample galaxies have Lb from 10 101 to 10 I(I9 L . 
On the other hand, Lism scatters from 10 39 4 erg/s to 10 424 erg/s. 

The right panel of Figure [3] shows the correlation between 
Lism within Ar e and Lg<x 2 , with <x denoting the central stel- 
lar velocity dispersion in each galaxy. If stellar motion is the 
main heat source for the hot ISM, its X-ray luminosity should 
be approximated by the input rate of the kinetic energy of the 
gas from stellar mass loss. This is proportional to L^cr 2 , be- 
cause the mass-lo ss rate is thought to be proportional to Lb (e.g., 
ICiotti et al.lll99ll) . Figure [3] also shows the expected energy in- 
put from stellar mass loss, La-, assuming a mass-loss rate of 



1.5 x 1O- U (L B /L ) f- 5 13 M G yr 1 dCiotti et all [l99ll) . Here, f 15 
is the stellar age in units of 15 Gyr, and we assumed a stellar age 
of 12 Gyr. Several galaxies have larger Lism than La-, while Lism 
of the other galaxies are similar to or smaller than L a . 

4.2.2. Temperature Profiles of the ISM and Classification 
with X-Ray Extended and X-Ray Compact galaxies 

Figure [4] shows the derived radial temperature profiles of the 
ISM. Some galaxies have gradually increasing temperature pro- 
files toward the outer radius. By contrast, other galaxies have flat 
or decreasing radial temperature profiles. The galaxies with pos- 
itive temperature gradients have high ISM temperatures of 0.8- 
1.5 keV at a radius of several times r e , which are comparable to 
those of galaxy groups. On the other hand, the temperatures of 
the other galaxies are systematically lower at 0.2-0.6 keV. The 
derived temperature profiles were fitted with the sum of a con- 
stant and single- or double-/? functions. 

The derived temperature pr ofiles are mostly c onsistent with 
previous results from ROSAT dMatsushital 1200 lb and Chandra 
( Fuk azawa et a n i2006trAthiv[ 12001 7 except central regions of 
several brightest galaxies with positive temperature gradients, 
due to higher angular resolution of Chandra. 

In Figure |5] we plotted kT(\-4r e )lkT{< lr e ) and kT(4- 
8r e )/kT(< lr e ) against L lSM /L B cr 2 . Here, kT(< lr e ), kT(l-4r e ), 
and £L(4-8r f ) correspond to the emission-weighted ISM tem- 
peratures of regions in r < \r e , \r e < r < 4r e and Ar e < r < 2>r e , 
respectively. In general, the temperature profiles of galaxies with 
Lism ^ La- have smaller kT{\-Ar e ) and kT{A-%r e ) than kT(< \r e ). 
By contrast, Lism of galaxies with kT(4-8r e )/kT(< \r e ) > 1.3 
are systematically larger than L tr . 

Our sample galaxies are divided into two types: X-ray faint 
galaxies with flat or negative temperature gradients and X-ray 
luminous galaxies with positive temperature gradients. The ISM 
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Fig. 2. (a) Innermost projected spectra observed with MOS (black) and PN (red), and innermost deprojected spectra with MOS 
(blue) and PN (magenta). These spectra are fitted with a vAPEC model plus power-law multiplied by the Galactic absorption (solid 
line). Dashed lines correspond to the contribution from each component. (b)(c) The actual data to model ratio from the fit in panel 
a. The meanings of colors correspond to that in panel a. (d) The same as panel b and c, but when using vAPEC + vAPEC model 
plus power-law multiplied by the Galactic absorption as fitting model. 
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Fig. 3. Lism(< 4r e ) of galaxies plotted against Lb (left panel) and Lgcr 2 (right panel). Symbols indicate galaxy categories defined in 
Section l4.2.2l for Xe galaxies (red crosses), field Xq galaxies (filled blue circles), and Xq galaxies in the clusters (open blue circles). 
The solid line represents the kinetic heating rate by stellar mass loss {La-). 
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Fig. 4. The derived temperature profiles of the ISM. The Xe galaxies are plotted in the top left panel, the Xc galaxies in clusters are 
in the top right, and the field Xc galaxies are in the bottom panels. Colors indicate individual galaxies. The solid lines represent the 
best-fit function. 



luminosities in the former type are consistent with heating by 
stellar motion, while galaxies of the latter type need additional 
sources of heating. Hereafter, we denote galaxies with kT(4- 
%r e )/kT(< \r e ) > 1.3 as X-ray extended (Xe) galaxies and oth- 
ers as X-ray compact (Xc) galaxies. The classification of each 
galaxy is summarized in Table Q] 



5. Spatial Analysis and Results 

5. 1 . X-ray surface brightness and gas density profiles 

We derived radial profiles of X-ray surface brightness from 
background-subtracted and vignetting-corrected X-ray images 
from MOS1 and MOS2. We considered only photons in the en- 
ergy band 0.8-2.0 keV, where the ISM emission dominates. PN 
data were not used for this analysis because of gaps between 
CCD chips. Figure [6] shows two representative X-ray surface 
brightness profiles, those of an Xe galaxy, NGC 4636, and an 
X c galaxy, NGC 720. 

Assuming circular symmetry, we deprojected the X-ray sur- 
face brightness profiles to derive gas density profiles. In order to 
subtract emission from outside the field of view, we first fitted 
the radial surface brightness profile within 8r e of each galaxy 
and assumed that the profile extends outside the field of view. 
The radial profiles of most of the galaxies were fitted with a /?- 
model. The brightness of several galaxies in clusters are clearly 
constant at the outer regions because of the surrounding intra- 
cluster medium (ICM). Therefore, we fitted these profiles with 
the sum of a /? model and a constant. Several X-ray luminous 
galaxies need a double-/? model to fit the surface brightness pro- 
files. Because at > 5r e , one of the /? models dominates, the pro- 
files at r > 5r e are fitted with a single-/? model. The derived 
density profiles are summarized in Figure|7] Since we need only 



the gradient of a density profile to derive a gravitational mass 
profile, the plotted density profiles were arbitrarily normalized. 

We also used Chandra data for 19 galaxies with sufficiently 
high signal-to-noise ratios to derive accurate X-ray surface 
brightness profiles within l-2r e . Radial profiles of X-ray surface 
brightness were derived from ACIS X-ray images in the energy 
band 0.3-2.0 keV. Then, we deprojected X-ray surface bright- 
ness profiles and derived gas density profiles in the same way 
as for the XMM-Newton data. As summarized in Figure [TJ nor- 
malized density profiles of the ISM derived from XMM-Newton 
and Chandra are mostly consistent with each other from Q.5r e to 
2r e . 

We then fitted the derived gas density profiles of XMM- 
Newton at r > 0.5r e and that of Chandra within r < 2r e of 
each galaxy simultaneously with a /? model, as f(R) = Sq(1 + 
(R/R C ) 2 )~ IN . Most of the density profiles were well fitted with 
this single-/? model (Figure |7). The derived R c values of galax- 
ies with Chandra data are almost about 0.05r e . Therefore, we 
fixed R c value to 0.05r e to fit the gas density profiles of galax- 
ies without Chandra data. For NGC 4552 we also fixed R c value 
to 0.05r e , bec ause of edge like struct ure at 0.4-0. 6r e observed 
with Chandra (Ma chacek et al.L [20061) . Several X-ray luminous 
galaxies need a double-/? model to fit the density profiles. 

The ISM and the hard component may have different surface 
brightness profiles. Therefore, we also derived gas density pro- 
files of galaxies whose Ljsm is smaller than 80% of the total lu- 
minosity within 4r e , directly from the normalization of the ISM 
component derived from spectral fittings. The best-fit /?-model of 
the density profile of each galaxy is plotted in Figure A. 1 . The 
best-fit values of /? derived in this way are mostly consistent with 
those derived from surface brightness profiles, although several 
galaxies show discrepancies of a few tens of %. 
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Fig. 5. The ISM temperature gradients parametrized by kT(l-4r e )/kT(< \r e ) (left panel) and kT(4-8r e )/kT(< \r e ) (right panel) 
plotted against Lism/^bct 2 - Meanings of the symbols are the same as those in Figure[3] The dotted line represents the rate of kinetic 
heating by stellar mass loss. The dashed lines indicate that kT{\-Ar e ) (or kT{A-%r e )) is 1 or 1.3 times higher than kT(< \r e ). 
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Fig. 6. X-ray surface brightness profiles derived from MOS images of NGC 720 (blue) and NGC 4636 (red). 




Fig. 7. ISM density profiles. The Xe galaxies are plotted in the top left panel, the Xq galaxies in clusters are in the top right, and the 
field Xq galaxies are in the bottom panels. These profiles are arbitrarily normalized. Colors indicate individual galaxies. Solid lines 
represent the best-fit function. Crosses and diamonds correspond to XMM-Newton and Chandra data, respectively. 
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Table 5. X-ray luminosities within 4r e of each sample galaxies 



Galaxy 


logLlSM" 


log L hmd h 


log L ISM a 




projection 


projection 


deprojection 




(erg/s) 


(erg/s) 


(erg/s) 


IC1459 


40.08 


40.28 


40.00 


NGC720 


40.46 


39.91 


40.38 


NGC1316 


40.63 


40.20 


40.55 


NGC1332 


40.12 


39.74 


40.04 


NGC1395 


40.43 


40.16 


40.34 


NGC1399 


41.13 


40.32 


40.98 


NGC1549 


39.39 


39.09 


— 


NGC3585 


39.36 


39.39 


— 


NGC3607 


40.37 


39.86 


— 


NGC3665 


40.36 


39.96 


— 


NGC3923 


40.74 


40.27 


40.68 


NGC4365 


39.58 


39.86 


— 


NGC4382 


40.25 


39.83 


— 


NGC4472 


41.40 


40.51 


41.36 


NGC4477 


39.96 


39.54 




NGC4526 


39.53 


39.64 




NGC4552 


40.62 


40.27 


40.60 


NGC4636 


41.46 


40.27 


41.44 


NGC4649 


41.00 


40.41 


40.99 


NGC5044 


42.49 


41.01 


42.41 


NGC5322 


40.14 


39.82 




NGC5846 


41.72 


40.56 


41.70 



"The X-ray luminosity of the thermal emission in the 

range of 0.3-2.0 keV derived from the projected 

and deprojected spectra. 
*The X-ray luminosity of the non thermal emission in 

the range of 0.3-2.0 keV derived from the projected 

spectra. 



6. Mass profiles 

We then calculated the total mass profile M(R) within a three- 
dimensional radius R from the obtained best-fit functions of ISM 
temperature T(R) and gas density n(R) profiles from the surface 
brightness, assuming hydrostatic equilibrium and circular sym- 
metry, by the equation 

JT(R).R ( dlnn(R) dlnT(R) \ 

Gum,, \ d\nR dlnR )' w 

where m p is the proton mass, k is the Boltzmann constant, G is 
the constant of gravity, and fi ~ 0.62 is the mean particle mass 
in units of m p . Figure [8] summarizes the derived mass profiles. 
For NGC 1549, NGC 4477, and NGC 5322, the mass profiles 
within 0.5r e were not plotted in the figure, since we used only 
XMM-Newton data for those galaxies. 

In addition, the upper and lower limits of the mass profiles 
were calculated considering the errors in the temperature and 
the temperature and density gradients of each data bin. The up- 
per and lower limits of the temperature and density gradients of 
the 2-th shell were obtained from the ratio of the value within 
i + 1-th shell to that within i - 1-th shell. Here, we used the 
temperature profiles derived from XMM-Newton. At r < \r e 
and r > lr e , density profiles derived from Chandra and XMM- 
Newton, respectively, were used. For NGC 4636, the density at 
0.3-0.4/v is significantly smaller than the best-fit function, due 
to the existence of complicated structure discovered by Chandra 
dJones et al.ll2002l) . Therefore, we ignored this shell when de- 
riving the mass profiles. As summarized in Figure [8] the total 
masses derived in this way are mostly consistent with those us- 
ing the best-fit functions. 



We also derived stellar mass profiles, using the deprojected 
de Vaucouleurs profile of iMellier & Mathe^ ( 119871) . assuming 
stellar Ml Lb is in the range from 3 to 8 in solar units. We plot- 
ted these profiles in Figure[8] Within \r e , the gravitational mass 
is consistent with the stellar mass. Further, the gradients of the 
gravitational mass are similar to those of the stellar mass. By 
contrast, outside the radius of a few r e , the derived gravitational 
mass becomes much larger than the stellar mass. These results 
indicate the existence of dark matter in the outer regions of early- 
type galaxies. 

When Lism is smaller than 80% of the total luminosity 
within 4r e , the total mass profile is calculated using the best-fit /3- 
model of the density profile from the spectral fitting. The results 
are compared with those from the surface brightness in Figure 
A.l. The two methods give similar total mass profiles within 
10-20%. For NGC 3585 and NGC 5322, these discrepancies are 
~30%, but within the large errors of the mass profiles of the two 
galaxies. Thus, the mass profiles of these galaxies obtained from 
two methods are consistent within the error. Hereafter, we use 
the total mass profiles derived from the surface brightness pro- 
files. 



7. Discussion 

The observed temperature profiles and X-ray luminosities of the 
ISM lead to a division of early-type galaxies into two categories: 
Xe galaxies and Xc galaxies. The Xe galaxies have increasing 
temperature profiles and Lism/^o- > 1, whereas the Xc galaxies 
have flat or negative temperature gradients and Lism/^o- ^ 1. 
Here, Lg- represents the expected energy input from stellar mass 
loss (iMatsushital |2001|) . In Section [TT] the derived ISM tem- 
peratures are compared with the stellar velocity dispersions. In 
Sections ^. 2l and l731 we derive gravitational-mass-to-light ratios 
in the B and K band, respectively, and constrain contributions 
from stellar mass and differences in dark mass between the Xc 
and Xe galaxies. Finally, in Section [7741 we discuss dark matter 
distribution in early-type galaxies themselves and their luminos- 
ity. 

7.1. ISM temperature vs. stellar velocity dispersion 

Figure|9]shows the correlation between central ISM temperature 
and central stellar velocity dispersion <x 2 (TableQTJi. The temper- 
ature roughly correlates with the stellar velocity dispersion. The 
parameter f3 spec denotes the ratio of stellar velocity dispersion 
to ISM temperature with /3 spec = fim p a 2 /kT, with fi indicating 
the mean molecular weight in terms of proton mass m p . For Xc 
galaxies, f} spec is about f} spec =0.5-1.0, which indicates that the 
ISM temperature s are c onsistent with heating due to stellar mo- 
tion ( Mat sushital 1200 lb . The galaxies with low cr tend to have 
low /3 values. It may be due to a selection effect, since brighter 
galaxies with higher ISM temperatures were the first proposed 
for observation. 

In Figure[10j we plotted the temperature profiles of the sam- 
ple galaxies scaled with central stellar velocity dispersion. There 
is no significant difference in the ISM temperature between the 
Xc and Xe galaxies at the central regions. Therefore, the ISM 
temperatures of the Xc and Xe galaxies would reflect the same 
potential. However, in the outer regions the Xe galaxies have 
higher ISM temperatures than the Xc galaxies for the same stel- 
lar velocity dispersion. This is due to the difference in potential 
due to the hot intra-group medium around the Xe galaxies. 
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Fig. 8. Integrated mass profiles of galaxies. The solid lines are total gravitational mass obtained from the best-fit functions of 
temperature and ISM density profiles. The dashed lines correspond to the upper and lower limits derived from the local gradients of 
temperature and density. We also plotted the stellar mass profiles assuming stellar M/L B to be 3 and 8 (dotted lines). 




2.3 2.4 
log a (km/s) 

Fig. 9. Central ISM temperatures plotted against central stellar 
velocity dispersion, <x. Colors indicate galaxy categories for Xe 
galaxies (red crosses), field Xc galaxies (filled blue symbols), 
and Xq galaxies in clusters (open blue symbols). For Xc galax- 
ies, the temperatures within 0.5 r e (circles) and \ r e (triangles) are 
plotted for the data with high and poor statistics, respectively. 
The dotted lines correspond to /5 spec of 0.25, 0.5, and 1.0 from 
the top to bottom. 



7.2. Mass-to-light ratio in B-band 

In order to compare the difference in dark matter profiles be- 
tween the Xc and Xe galaxies, we used the mass-to-light ratio, 
Ml Lb- We assumed the de Vaucouleurs law for the stellar dis- 
tribution. Figure [TT1 summarizes the integrated M/Lb profiles. 
For NGC 1549, NGC 4477, and NGC 5322, profiles within 0.5r e 
were not plotted in the figure, since we used only XMM-Newton 




1 

R/r. 



Fig. 10. Temperature profiles of the ISM scaled with central stel- 
lar velocity dispersion. The solid red, solid blue, and dashed blue 
lines represent the Xe galaxies, field Xc galaxies, and Xc galax- 
ies in clusters, respectively. 



data for those galaxies. From 0.2r e to lr e , M/Lb of each Xc 
galaxy is nearly constant at 3-10, which is consistent with typi- 
cal values of stellar M/Lb- Even in the Xe galaxies, M/Lb within 
0.5-1 r e are also flat. These results suggest that stellar mass dom- 
inates within r < 0.5- lr e . In NGC 4636, a wavy M/Lb profile is 
seen at r < Q.3r e , which is artificially caused by the complicated 
X-ray structures i n the central regions discovered by Chandra 
dJones et al.l2002l) . Hereafter, we denote integrated M/Lb within 
sr e as M/Ls(r < sr e ), where s is a constant. 

The top left panel of Figure [12] shows the relationship 
of M/Ls(r < Q.5r e ) to total B-band luminosity, Lb- In this 
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Fig. 8. (continued) 



plot, we add M/Lb of the cD galaxy of the Virgo cluster, M 
87 ob tained from XMM-Newton observation dMatsushita et all 
120021) . M/L B (r < 0.5r e ) of the X E and X c galaxies is about 
10, with significant scatter. A correlation between M/Lb and 
Lb elliptical galaxie s was found by optical measurements 
dGerhard et all 120011) . where gravitational mass was derived 
from stellar velocity dispersion at the central region of elliptical 
galaxies. Our M/Lb and L B re lation scatters around the relation 
found bv lGerhard etaH d200il) . 

On the other hand, at r > lr e , the derived M/Lb starts to in- 
crease (FigurefTTT). The Xc galaxies have similarly shaped M/Lb 
profiles. M/L B (r < 3r e ) and M/L B (r < 6r e ) of the X c galax- 
ies are about 6~ 25M /L Q . The Xe galaxies have systematically 
larger M/L B values than the Xc galaxies at > 3r e . (Figure [12b. 
The Xe galaxies have M/L B (r < 3r e ) and M/L B {r < 6r e ) of 
25 ~ 4OM /L G and 40 ~ 100M o /L o , respectively. These re- 
sults indicate that dark matter is common in early-type galaxies. 
In addition, the X E galaxies have more dark matter than the Xc 
galaxies. 



7.3. Mass-to-light ratios in Kband 

Historically, the B-band mass-to-light ratio has been used in 
such studies. However, K-band luminosity well describes the 
stellar mass. Thus, we also derived the K-band mass-to-light 
ratio to study the dark matter profiles. We calculated the K- 
band luminosity Lk from the Two Micron All Sky Survey 
(2MASS). The effect of Galactic extinction was corrected us- 
ing the NASA/IPAC Extragalactic Database (NED). The rela- 
tionship between L B and L K is close to that of B - K = 4.2, 
which is the appro priate color of stars in early-type galaxies 
dLin&Mohil 12004 . with some scatter (Figure [13V 



The right panel of Figure QT| shows the integrated profiles of 
the K-band mass-to-light ratio, M/Lk, of the sample galaxies. At 
the central region, the Xe and Xc galaxies have similar M/Lk at 
~ 1. 

In Figure [14] we compared the relationship of M/Lx(r < 
0.5r e ), M/L K (r < 3r e ), and M/L K (r < 6r e ) with L K except M 
87. For early-type galaxies, the stellar K-band mass-to-light ra- 
tio, M/L K , is about 1M /L Q . The scatter of M/L K (r < 0.5r e ) 
values becomes smaller than that of M/Ls{r < Q.5r e ) values 
(Figure [T4V Both the Xe and Xc galaxies have M/Lk values of 
~ 1M O /L , and no correlation with Lk- This result indicates that 
stars dominate the mass within Q.5r e , and we observed the stel- 
lar M/Lk in these region, except the cD galaxy M 87. M 87 may 
contain similar amount of dark matter with stellar mass within 
0.5r e , since M/Lk is a factor of 2 larger than those of Xe and 
Xc galaxies. On the other hand, the M/L^ir < 3r e ) values of 
the Xe and Xc galaxies are 4 ~ 10 and 1 ~ 4M e /L Q , respec- 
tively (Figure [l4j. The M/L K (r < 6r e ) values of the X E and X c 
galaxies are ~ 10 and 2 ~ 6M O /L , respectively (Figure [T4V 

7.4. Dark matter in early-type galaxies 

We normalized the M/L K profile with M/L K (r < 0.5r e ) for 
each galaxy (Figure |T5T>. Then, the profiles of mass-to-light ra- 
tios of Xc galaxies have similarities. We also derived the ratio 
of M/L K (r < 3r e ) and M/L K (r < 6r e ) to M/L K (r < 0.5r e ) for 
each galaxy. The ratios of M/L K {r < 3r e ) and M/L K {r < 6r e ) 
to M/Lx(r < Q.5r e ) are similar among the Xc galaxies at 2 and 
3 ~ 4, respectively (Figure [T6T >. Considering that stellar mass 
dominates within \r e , and assuming that the stellar M/Lk is 
nearly constant within a galaxy, the Xc galaxies contain simi- 
lar amounts of dark matter. Within 3r e , the mass of dark matter 
is similar to the stellar mass, while within 6r e , the former is 2-3 
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Fig. 11. Profiles of integrated M/Lg (left) and M/Lk (right). The meanings of colors and lines are the same as those in Figure [TOl 
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Fig. 12. M/Ls(r < 0.5r e ) (top left), M/Ls(r < 3r c ) (top right), and M/Ls(r < 6r e ) (bottom) against Lb- The quadrangle represents 
the value of M87 bv lMatsushita et all ([2002). Meanings of ot her symbols are the sa me as those in Figure[3] The solid line represents 
the correlation derived from the stellar velocity dispersion bv lGerhard et alj d200ll) . 



times larger than the latter. These ratios should reflect the poten- 
tial of early-type galaxies themselves. Thus, the dark mass distri- 
bution in early-type galaxies is slightly more extended than that 
of stars. It is thought that early-type galaxies have 10 times more 
dark m ass than stellar mass as in spiral galaxies (e.g JCiotti et alj 
H99l . However, the galaxies themselves may not contain such 
large amounts of mass, at least within several times r e . This re- 
sult should be important for the study of the origin of the dark 
matter content in early-type galaxies and for the study of the for- 
mation and evolution of these galaxies. 

By contrast, the Xe galaxies have systematically larger ra- 
tios; i.e., they have more dark matter in their outer regions. 
Figure[l2]shows that the galaxies with a larger LxlL a have more 
dark matter at r > lr e . The ratio of M87 may also in the rela- 
tion between Lx/L a and dark to stellar mass ratio, considering 
that M 87 may contain similar amount of dark mass to stellar 
mass within Q.5r e . In other words, the X-ray luminosity of the Xe 
galaxies may be determined in relation to their potential struc- 



ture, as indicated by iMatsushital d2001l) and iMatsushita et al.l 
(120021) . The difference in temperature profiles between the Xe 
and Xq galaxies would be due to differences in the surround- 
ing gravitational potential. These results suggest that the X-ray 
luminous early-type galaxies commonly sit in the center of a 
large-scale (a few hundred kpc) potential well, which leads to 
their high luminosities. Other galaxies may lack such a large- 
scale potential well and contain their own dark matter. 

8. Conclusion 

We analyzed 22 early-type galaxies using XMM-Newton and 
Chandra data. To derive the gravitational mass profiles, we ob- 
tained the temperature and ISM density profiles through spec- 
tral fitting and spatial analysis, respectively. We classified the 
galaxies into two categories, Xe and Xq galaxies, on the basis of 
whether the temperature gradient is positive or negative toward 
the outer radius. The ISM luminosity of the Xq galaxies is con- 
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Fig. 14. Integrated M/L K (r < 0.5r e ) (top left), M/L K (r < 3r e ) (top right), and M/L K (r < 6r e ) (bottom) against L K . Meanings of the 
symbols are the same as those in Figure [T2l 
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Fig. 16. Ratio of the integrated M/L K (r < 3r e ) and M/L K (r < 6r e ) to M/L K (r < 0.5r e ) plotted against L K . Meanings of the symbols 
are the same as those in Figure [T2l 
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Fig. 17. Ratio of the integrated M/L K (r < 3r e ) and M/L K (r < 6r e ) to M/L K (r < 0.5r e ) plotted against Lism/L b ct 2 . Meanings of the 
symbols are the same as those in Figure[T2] The dotted line represents the kinetic heating rate by stellar mass loss. 
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Galaxy 


M/L B (< Q.5r e ) 


M/L B (< 3r e ) 


M/L B (< 6r c ) 


M/L K (< 0.5r e ) 


M/L K (< 3r e ) 


M/L K (< 6r e ) 




(MJLe) 


(M /L ) 


(M G /L G ) 


(M G /L G ) 


(M G /L G ) 


(M G /L G ) 


IC1459 


8.8 


18.1 


29.9 


1.4 


2.9 


4.8 


NGC720 


8.7 


19.5 


33.6 


1.7 


3.9 


6.7 


NGC1316 


9.1 


17.8 


29.5 


1.5 


3.0 


4.9 


NGC1332 


10.6 


16.5 


27.3 


1.7 


2.7 


4.4 


NGC1395 


10.9 


27.7 


50.9 


1.5 


3.9 


7.2 


NGC1399 


14.3 


27.1 


58.1 


1.6 


3.0 


6.4 


NGC1549 


10.2 


20.9 


— 


1.7 


3.5 


— 


NGC3585 


4.0 


6.9 


11.4 


0.7 


1.2 


1.9 


NGC3607 


8.1 


17.3 


30.8 


1.5 


3.2 


5.8 


NGC3665 


4.8 


10.0 


16.6 


0.9 


1.8 


3.0 


NGC3923 


8.2 


13.7 


20.4 


1.1 


1.8 


2.6 


NGC4365 


8.4 


17.4 


— 


1.6 


3.2 


— 


NGC4382 


3.3 


6.0 


8.7 


0.7 


1.2 


1.8 


NGC4472 


6.5 


23.8 


40.4 


1.3 


4.7 


7.9 


NGC4477 


10.0 


22.0 


31.0 


2.0 


4.3 


6.1 


NGC4526 


5.3 


9.2 




0.9 


1.5 




NGC4552 


7.4 


10.9 




1.3 


1.9 




NGC4636 


9.8 


43.0 


75.3 


1.7 


7.3 


12.8 


NGC4649 


9.1 


21.4 


36.4 


1.6 


3.8 


6.5 


NGC5044 


6.9 


41.5 


109.9 
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Fig. 13. Relationship between B-band and K-band luminosity of 
the sample galaxies. Meanings of the symbols are the same as 
those in Figure[3] The dotted line correspo nds to the appropria te 
color B - K - 4.2 for early-type galaxies dLin & MohrM2004l) . 



sistent with the energy input from stellar mass loss. By contrast, 
the Xe galaxies have larger ISM luminosity. 

At the central regions, r < 0.5 - \r e , the derived integrated 
M/Lk of both of the Xq and Xe galaxies are about 1 and have 
smaller scatter than M/Lb- The values and profiles of M/Lk indi- 
cate that stellar mass dominates the total mass in these regions. 
In the outer regions, M/Lb and M/Lk of the Xe galaxies are 
higher than those of the Xc galaxies. 

On the basis of these results, we can conclude the follow- 
ing. The normal early-type galaxies, Xq galaxies, contain their 
own dark matter at amounts that are 1.5~2 times larger than the 
stellar mass within 5r e , The Xe galaxies are located as the cen- 
tral galaxy in a larger scale potential structure, such as a galaxy 
group. This fact directly affects the gravitational potential profile 




Fig. 15. Profiles of M/L K normalized b y M/L K (r < 0.5r ? ). We 
also plotted the M/L K profile of M87 bv lMatsushita et al.l(l2002l) 
(orange solid line). The meanings of other colors and lines are 
the same as those in Figure [TOl 



of the galaxy itself, and causes it to contain significantly higher 
amounts of dark matter than that in the Xq galaxies. This differ- 
ence in the gravitational potential leads to the difference in the 
temperature profile and X-ray ISM luminosity between the Xe 
and Xq galaxies. 
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Galaxy ring" kT b O Si Fe logLjsM'' log L/ Ulr / J ^--/d.o.f. 

(r e ) (keV) (solar) (solar) (solar) (erg/s) (erg/s) 

IC1459 51)75 0.606+™;; 0.56+°;* 0.54+° 44 0.36+°°= 39760 4O07 242 / 201 

0.5-1 0.622+g : g|| — ' — " — ' 39.32 39.54 — 

1- 2 0.574+™P — — — 39.47 39.42 — 

2- 4 0.589+g : g|| — — — 39.48 39.12 — 

4-8 0.638+™g — — — 39.68 38.99 — 

NGC720 51 - 601 S 0.49+™^ 0.3 1+°:; 4 0.50+°°= 40761 3946 207/ 157 

1- 2 °- 564 -8np 39-93 39.23 

2- 4 0.554+™'' — — — 40.01 39.56 — 

4-8 0.480+™g — — — 39.94 39.54 — 

NGC1316 OU25 0.728+™™ 0.48+™= 0.24+°° 4 0.36+°°j 4O09 39745 1602/939 

0.25-0.5 0.632+™ 07 , 39.82 39.31 

0.5-1 0.599+™°;? 40.00 39.61 

1- 2 0.664+™™ — 39.87 39.35 

2- 4 0.373+™]= 39.81 39.68 

4-8 0.616+™^ — — — 39.90 39.77 — 

NGC1332 OHl 0.6 17+™ 0.72+™^ 0.89+°]° 0.56+™= 39~72 39714 337 / 227 

0.5-1 0.580+™]' — — ' — ' 39.36 39.06 — 

1- 2 0.486+™'= — — — 39.40 39.15 — 

2- 4 0.495+™™ 39.49 39.18 

4-8 0.542+™^ — — — 39.34 38.89 — 

NGC1395 (M)l 0.597+™™ 0.79+°,;= 1.15+J^ 0.80+°°^ 39763 39734 381 /342 

0.5-1 0.607+™" — ' — ' — ' 39.67 39.37 — 

1- 2 0.710+™!° — — — 39.89 39.69 — 

2- 4 0.748+™'° — — — 40.02 39.71 — 

4-8 0.762+™ — — — 40.02 39.94 — 

NGC1399 (M)l 57846+™ 0.72+",j 4 1.24+°|= 1.07+°™ 4052 39756 374 / 250 

0.5-1 0.987+™]= — " — ' — ' 40.36 39.44 — 

1- 2 L169 -8p 40 - 41 3993 

2- 4 1.342+g°TT 40.73 39.77 

4-6 1.440+™]° — — — 40.62 39.85 — 

6-8 1.535+™! ~ — 40.43 39.87 — 

NGC1549 51 0.398+™™ 6725+™ lW' 0.27+°°° 3911 39l7 93/ 111 

1-4 0,441+°$! — ' — — ' 39.03 37.73 — 

NGC3585 51 0.412+J™ 0.78+™° fTOCT 1.03+jJ-^ 1 38769 7397T2 241 / 159 

1- 2 0.382+° : J , > — ' — — ' 38.70 39.00 — 

2- 4 0.338+||i — — — 39.11 38.14 — 

4-8 0.399+™g — — — 39.41 39.06 — 

NGC3607 W 0.595+°| 4 0.09+™^ 0.18+°* 0.24+J™ 39768 39751 177/ 193 

1- 2 0.550+™| — ' — ' — ' 39.73 39.26 — 

2- 4 0.566+™g — — — 40.13 39.34 — 
4-8 0.449+° : °J= — — — 40.20 40.18 — 

NGC3665 (HI 0.448+™* 0.15+™= 0.01+° 4 | 0.23+°°^ 39747 39741 275 / 223 

0.5-1 0.382+™;= — ' — ' — ' 39.65 39.29 — 

1- 2 0.378+™;; — — — 39.82 39.38 — 

2- 4 0.377+™;; — — — 39.94 39.34 — 

4-8 0.435+™" — — — 39.82 O00 — 

NGC3923 61)725 0.613+°™ 0.37+™= 0.24+°]J 0.37+°°; 40~T6 397l2 319/220 

0.25-0.5 0.559+™ — ' — ' — ' 40.03 39.27 — 

0.5-1 0.479+™]; — — — 40.04 39.52 — 

1- 2 0.491+™' — — — 40.00 39.65 — 

2- 4 0.562+™'' — — — 39.91 39.88 — 

4-8 0.471+° : ]^ — — — 39.77 39.83 — 

NGC4365 51 0.522+™| 1.00+°™ 0.96+^ 4 6747+JTTJ WE 39747 177/ 190 

1-4 0.624+°™ — ' — ' — ' 39.37 39.64 — 

NGC4382 51X5 0.404+™°= 0.17+™^ 57T6+™ 0.33+°°° 39~36 39720 244/211 

0.5-1 0.384+™]= — ' — ' — ' 39.53 39.09 — 

1- 2 0.349+™]] — — — 39.76 39.07 — 

2- 4 0.315+° : ™ — — — 39.80 39.45 — 
4-8 0.266+™=' — — — 39.70 39.43 — 

NGC4472 0-0.25 0.782 5777 TTB 6792 45761 39753 782/272 

0.25-0.5 0.903 0.75 0.97 0.77 40.47 39.49 964/216 

0.5-0.75 0.995 1.01 1.53 1.11 40.40 39.58 667/229 

0.75-1 1.031 0.86 1.68 1.10 40.36 39.52 541 / 195 

1- 1.5 1.044 0.94 1.76 1.22 40.59 39.76 808/ 194 
1.5-2 1.160 0.58 1.38 1.06 40.50 39.65 625/ 186 

2- 3 1.227 0.58 1.35 1.00 40.61 39.78 427/ 193 

3-4 1.270+™]; °j3g!n!| i- 27 !"?" °- 67 -w 4ft38 3942 201 / 115 

NGC4477 51 67446+™= 0.25+™;! LOW 5717™* 39771 39761 97 / 79 

1- 2 0.452+™ 4 39.56 38.60 

2- 4 0.273+° : ™ — — — 38.67 39.32 — 

4-8 0.122+° : g] — — — 38.99 39.28 — 

NGC4526 51)1 0.353+™^ 0.29+™* LOlF 0.40+°°° 39721 39729 77/117 

0.5-1 0.322+™| — ' — — 38 ' 92 3892 — 

1-4 0.280+°™ — — — 38.97 39.20 — 

NGC4552 51X5 0.679+™;° 0.52+™= 0.56+° |= 0.41+°°^ 45721 39781 275/ 182 

0.5-1 0.582+™,!,° 40.08 39.54 

1- 2 0.500+™]= — — — 39.97 39.63 — 

2- 4 0.528+™;= — — — 39.66 39.63 — 
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13 +0j * 
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"The inner and outer radii to integrated the spectrum. 
h kT is ISM temperature. Lism and Lj lcirc i are the X-ray luminosities of the thermal emission and the non 

thermal emission in the range of 0.3-2.0 keV, respectively. 
* fixed to the value. 
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Table 4. Spectral fitting results of the deprojected annular spectra. 
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"The inner and outer radii to integrated the spectrum. 

b kT is ISM temperature. Ljsm ar| d are me X-ray luminosities of the thermal emission and the non 

thermal emission in the range of 0.3-2.0 keV, respectively. 
c We fitted the spectra with 2T model. 
* fixed to the value. 
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«/'. «/r. R/r. 

Fig. A.l. (Left) The ISM temperature profiles of sample galaxies. The cross data points represent temperature profile derived from projected 
spectra, and the diamond data points are from deprojected spectra. The solid lines represent each best-fit function. The black data points are used 
to fit the function, but the red data points are not. (Center) The ISM density profiles. The cross data points represent XMM-Newton data, and the 
diamond data points are Chandra data obtained from spatial analysis. The black solid lines represent each best-fit function. Meanings of the red 
and black colors are the same as left panel. The blue solid lines show best-fit function of ISM density profile obtained from spectral fit. These 
ISM density profiles are normalized by the value at \r e . (Right) The integrated mass profiles. The solid lines are total mass obtained from the 
best-fit function of temperature and ISM density profiles from the surface brightness (black) and spectral fit (blue). The dashed lines are upper 
and lower limit. We also plotted the stellar mass profiles assuming stellar M/L B to be 3 and 8 (red dotted lines). 
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Fig. A.l. (continued) 
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Fig. A.l. (continued) 
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Fig. A.l. (continued) 



